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Abstract. 

In suitably chosen domains of space-time, the world function may be a 
powerful tool for modelling the deflection of light and the time/frequency transfer. 
In this paper we work out a recursive procedure for expanding the world function 
into a perturbative series of ascending powers of the Newtonian gravitational 
constant G. We show rigorously that each perturbation term is given by a line 
integral taken along the unperturbed geodesic between two points. Once the world 
function is known, it becomes possible to determine the time transfer functions 
giving the propagation time of a photon between its emission and its reception. 
We establish that the direction of a light ray as measured in the 3-space relative 
to an observer can be derived from these time transfer functions, even if the 
metric is not stationary. We show how to derive these functions up to any given 
order in G from the perturbative expansion of the world function. To illustrate 
the method, we carry out the calculation of the world function and of the time 
transfer function outside a static, spherically symmetric body up to the order G^ , 
the metric containing three arbitrary parameters /3, 7, S. 



PACS numbers: 04.20.Cv, 04.25.-g, 04.80.-y 



§ To whom correspondence should be addressed (leponcin@danof.obspm.fr) 



World function and time transfer 



2 



1. Introduction 

With advances in technology, it wih become indispensable to determine some of the 
relativistic effects in the propagation of light beyond the first order in the Newtonian 
gravitational constant G, in particular in the area of space astromctry. The Global 
Astrometric Interferometer for Astrophysics (GAIA) [1] and the Space Interferometric 
Mission (SIM) [2] are already planned to measure the positions and/or the parallaxes 
of celestial objects with typical uncertainties in the range 1-10 /iarcsccond (/ias). The 
Laser Astrometric Test Of Relativity (LATOR) mission is an even more ambitious 
program since its primary objective is to measure the bending of light near the Sun 
to an accuracy of 0.02 /las [3]. With such requirements in experimental precision, it is 
clear that the effects of the second order in G must be taken into account in the near 
future since they result in several /zas deflection for a light ray grazing the Sun [4-7] . 

In order to obtain a convenient modelling of the above-mentioned projects, it 
is necesssary to determine the deflection of light rays between two points xa and 
xb of space-time. In almost all of the theoretical studies devoted to this problem, 
the properties of light rays arc determined by integrating the differential equations 
of the null geodesies. This procedure is workable as long as one contents oneself 
with analyzing the effects of first order in G, as it is proven by the generality of the 
results obtained in [8-15]. Unfortunately, analytical solution of the geodesic equations 
requires cumbersome calculations when terms of second order in G are taken into 
account, even in the case of a static, spherically symmetric space-time [4-6, 16]. The 
goal of the present paper is to develop a procedure allowing to determine the deflection 
of light and the travel time of photons without integrating the geodesic equations. 

We assume that space-time is globally regular with the topology M x M^, i.e. 
without event horizon. According to these assumptions, we admit that the Lorentzian 
metric is represented at any point x by a series in ascending powers of G 

oo 

9^Ax,G)=9j°J + Y,G''9^;i\x), (1) 

71=1 

where 

alS = Vtiu = diag (1, -1, -1, -1) 

in any quasi Cartesian coordinates. Neglecting all terms involving g\xl^^\ g]!i^'^^ ■, ■■■ 
defines the so-called nth post-Minkowskian approximation. 

Essentially, our method uses the notion of world function as developed by Ruse 
[17] and Synge [18, 19]. The world function fl{xA,XB) is defined as half the squared 
geodesic distance between two points xa and x b ■ Before entering into the presentation 
of the method, it must be noted that this definition raises a serious problem, since 
there could exist points xa and xb connected by several geodesies. Given a point xa, 
it is usual to define the normal convex neighbourhood of xa, say Nxai the set of 
all points x joined to xa by a unique geodesic path [20]. The world function Q.{xa, x) 
is then well defined for x G J^xa- agreement with the assumptions mentioned in 
the previous paragraph, we suppose that if xb is joined to xa by a null geodesic, 
then this geodesic is unique. Denoting the light cone at xa by Cx^^ this means that 
CxA C NxA ■ Thus we can use the world function to study the deflection of the light and 
the travel time of photons between xa and x b ■ Let us point out that our assumptions 
are compatible with the existence of multiple-images effects involving different travel 
times for the photons (see, e.g., [21] and references therein). 
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We shall see in section 2 that the invariant two-point function fl(xA, xb) has 
very interesting properties. First of all, considering the geodesic path connecting xa 
and xb, it may be shown that the tangent vectors at xa and xb are provided by the 
gradients of the world function. Second, it is clear that xa and xb are joined by a null 
geodesic if and only if ^}{xa,xb) = 0. As a consequence, the deflection of a light ray 
emitted at a finite distance may be easily determined when the world function of a 
given space-time is known, just as the travel time of a photon between two points. This 
means that the world function allows to calculate the time/frequency transfer functions 
and can be extremely useful in situations involving a comparison between two distant 
clocks, as it has been shown by two of us in a recent work [22] extending the results 
obtained in [23] and [24]. Other possible applications deserve to be emphasized. For 
example, the quantum phase accumulated by a freely falling particle is proportional 
to the geodesic distance covered through space-time [25]. So, the world function can 
be relevant for modelling situations like neutrino oscillations [26]. 

Taking equation Q into account, it is possible to assume that Q{xa, xb) admits 
an expansion as follows 

oo 

nixA,XB) = n<'°\xA,XB) + J2G^n<-^^ {xa,Xb)- (2) 

71=1 

The primary objective of the paper is to yield a recursive procedure for determining 
each term n'--'^\xA, xb) within the nth post-Minkowskian approximation without 
calculating the geodesic joining xa and xb- 

Our method is inspired by a procedure found by Buchdahl [27] for obtaining 
the perturbative expansion of the characteristic function defined a,s V = -y/2 \ 
However, the gradients of V{xa, xb) are not defined when xa and xb are linked by a 
null geodesic. So we prefer developing a new procedure involving Q itself. Considering 
that XA and xb are fixed points, we use the Hamilton- Jacobi equation satisfied by 
fl{xA, x) at point x. It is easily seen that this equation may be replaced by an infinite 
set of ordinary differential equations for the perturbation terms il^"'\xA, x) when x 
is constrained to move along the zeroth-order geodesic connecting xa and xb- As a 
consequence, each term il^"^{xA,XB) is obtained in the form of a line integral along 
a straight line in the background Minkowski metric. Thus, we obtain a recursive 
procedure which completely avoids the calculation of the perturbation of the geodesic 
joining the given end points. 

Within the first post-Minkowskian approximation, we recover a simple integral 
expression of Jl*-^-* already given by Synge [28]. In contrast, the fact that one 
has to calculate line integrals with integrands containing products of the first-order 
derivatives of terms n^^^~P\ p = 1, n — 1 proves to be awkward in attempts at 
getting explicit expressions of fi^") for n > 2. Fortunately, at least in the case where 
n = 2, we obtain a simplified formula for fJ*^^) where the quadratic terms in the first 
derivatives of fl^^^ appear only as boundary terms. Then, we develop an iterative 
procedure allowing to deduce the propagation time of a photon correct to within 
(n + l)th-ordcr terms from the expressions 

Finally, we derive the world function and the propagation time of light in a static, 
spherically symmetric space-time within the second post-Minkowskian approximation. 
In order to get expressions valid for a wide class of gravitational theories, we 
assume that the metric involves three arbitrary parameters: the usual post-Newtonian 
parameters f3 and 7 [29], and a supplementary parameter 5, introduced to take into 
account the term of order in spatial components gij (see, e.g., [30] and references 
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therein). Complete results were obtained by John [31] involving quite cumbersome 
calculations. Without any reference to this earlier paper, a complicated integration 
of the geodesic equations was performed by Richter and Matzner [16] in order to 
calculate the time delay for a special coordinate system adapted to the direction 
of the light ray at the emission point. More recently, Brumbcrg [32] has used the 
Euler-Lagrange equations of null geodesies to find the propagation time of light as a 
function of the boundary conditions in the second post-Minkowskian approximation of 
the Schwarzschild metric. For our part, we do not need determining the null geodesies 
to get the world function and the propagation time of light and we use only elementary 
quadratures to perform the calculations. 

The paper is organized as follows. In section 2 the general definition of the world 
function and its fundamental properties are recalled. In section 3 we introduce two 
distinct time transfer functions and we give their general relations with the world 
function. We establish that all the calculations relative to the direction of a light 
ray or to the frequency transfers can be performed when one of the two time transfer 
functions is explicitly known. In section 4 we set out the recursive method which 
allows to determine the world function in the nth post-Minkowskian approximation 
for any n > 2. In section 5 we show how to deduce the time transfer functions from 
the world function at any post-Minkowskian order of approximation. In section 6 
we determine explicitly the world function and the time transfer outside a static, 
spherically symmetric body within the second post-Minkowskian approximation. 
Some concluding remarks are given in section 7. In Appendix A we give an alternative 
derivation of the integral expression for fl^'^\xA,XB), based on an integration of the 
Euler-Lagrange equations of geodesies. 

In this paper c is the speed of light in a vacuum. The Lorentzian metric of space- 
time V4 is denoted by g. The signature adopted for g is (H ). We suppose that 

space-time is covered by some global quasi-Galilean coordinate system {x'^) = (a;°, x), 
where = ct, t being a time coordinate, and x = (x'). We assume that the curves 
of equations .x* = const are timelikc, which means that goo > anywhere. We employ 
the vector notation a in order to denote {a^,a'^,a'^) = (a*). Considering two such 
quantities a and b we use a • 6 to denote o'6' (Einstein convention on repeated indices 
is used). The quantity |a| stands for the ordinary Euclidean norm of a. For any 
quantity f{x^), f^a denotes the partial derivative of / with respect to x". The indices 
in parentheses characterize the order of perturbation. They are set up or down, 
depending on the convenience. 

2. Definition and fundamental properties of Synge's world function 

Since Synge's world function is not widely employed, we briefly recall how its 
fundamental properties can be straightforwardly derived from the variational principle 
defining the geodesic curves (for the characteristic fimction V = ^2 | $1 1, see e.g. 
Buchdahl [33-35]). We begin with recalling some useful results concerning the general 
variation of a functional defined as an action in classical mechanics. 

Let XA and a:;^ be two points of space-time connected by a differentiable curve 
Cab defined by the parametric equations x°' = x°'{Q, with (a < C < Cb- Given a 
Lagrangian function L{x",x^) with x^{Q = dx^{Q/d(, one can define the functional 
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S[Cab]= r L{x^{(:lxP{0)dQ, (3) 

the integral being taken along Cab- Let us consider now an arbitrarily neighbouring 
curve C'j^ig, connecting points xa' and xb', represented by parametric equations 
x'" = a;'"(C'), with Ca' < C < Cb'- Let x'X + Sx^ = x% and x% + Sx% = x%, 
be the coordinates of xa' and xb' respectively. Putting Sx"{() = x'"(C) — x"{Q, one 
can define a one-to-one correspondance between Cab and Ca'B' by the infinitesimal 
transformation 

x'^iC) ^ x"^iC + SO = x^iO + fe"(C) , (4) 

where 

5x"(C)=fe"(C)+i"(C)<5C (5) 

with the boundary conditions 

Sx''{Ca)^Sx'X, 5x''{Cb)^Sx%. (6) 

Performing an integration by parts leads directly to the following expression for 
the quantity SS — S[Ca'B'] — S[Cab] (see, e.g., [36]) 



Ca 



'dx^ ~ dC \di^ 



Sx'^iOdC+iPaSx'^ ~HSC]a , (7) 



where pa is the 4-momentum belonging to x°' and H is the Hamiltonian: 
dL 

^"^ai^' H=p^i"-L. (8) 

In the following we assume that there exist domains V of space-time such that 
whatever xa (zV and xg (zV, xa and x^ are linked by a unique curve Tab realizing 
an extremum of the functional S. Then, describing Tab by parametric equations 
x" = x"{C)i we can associate to any domain V a function S of xa, xb, Ca, Cb defined 
as 

S{xa,xb,Ca,Cb)= f L{x"{C),i^{C))dC. (9) 

Extending a terminology currently used in mechanics, we shall call S the characteristic 
(or principal) function belonging to the Lagrangian L adapted to the domain T). Since 
the functions x" (C) satisfy the Euler-Lagrange equations 

dl^_d_ / dL 
dx°' ~ dC \di^ 

the total variation of S{xa,xb, Ca, Cb) reduces to the boundary terms in equation 0. 
Consequently, one has the relations 

{Pa)A , TTTT = (Po')B , (H) 



0, (10) 



dx'^ ' dxg 
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Considering that _ff is a function of x" and p^, then substituting for {pa)A E^nd {pa)B 
from equations ((TT|) into equations ifT^ . it is easily seen that S{xatXb,C,a,C,b) satisfies 
a Hamilton- Jacobi equation both at x'^ and at x'^ 

Let us apply these results to space-time V4 endowed with a Lorentzian metric 
and insert in equation Q the Lagrangian defined by 

L = ^9a0i"x>' . (14) 

Each curve Tab parametrized by C for which S has an extremum is a geodesic path 
joining xa and xb- Parameter C, is said to be afhne. According to |(Hl, we have now 

Pa = 9al3x'^ , H = ]^g°'^PaPp , (15) 

from which it is easily deduced that H = L on any geodesic path. Owing to the fact 
that L does not contain C, as an explicit variable, this last equality implies that the 
Lagrangian L defined by (|14|l is a constant of the motion. 

We are only concerned here with a weak-field metric represented by equation 
throughout space-time. So we henceforth restrict our attention to the domains 
= AfxA^ being defined in Introduction. Then whatever xb G ^xai there 

exists a unique geodesic path V ab connecting xa and xb- Let us denote by the 
characteristic function belonging to the Lagrangian (|14|l adapted to Mxa ■ According 
to (O we have 

?L(xa,xbXaXb) = \ gaf3x'^i'^dC, (16) 

^ Jca 

the integral being taken along Tab- Since L ^ H is a constant of the motion, we have 
L ~ Ha = Hb on T ab- Inserting these relations into equation H16(l yields 

TT TT ^{xa,xb,Ca,Cb) 

HA = tiB = . (17) 

Qb - Qa 

Substituting these expressions of Ha and Hb into equations ifT^ , we get 
dVt _ Vl dVL _ Q. 

dU^ Cb-Ca' dC^~~CB-W ^ ^ 

Integrating these equations, we find that may be written as 

r>i /- /- \ ^{xa,xb) , . 

il[XA,XB,CAXB) ^ —: T — , (19) 

Cb - Ca 

where fl{xA,XB) is a function of xa and of xb- This two-point function, which is 
symmetric in xa and xs, is called the world function. Denoting by A the unique affine 
parameter such that Xa = and As = 1, we infer from equations (|16|l and IjlQI) that 
fl{xA, Xb) may be written as 



1 dx^^ dx^ 

^{xa,xb) = - 3^iAx'"{X)) — —d\, 



(20) 



the integral being taken along Tab- It can be easily seen that 



VL{xa, Xb) = ^eAB[sABf , (21) 
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where cab — l;Oj~l when T^b is a timehke, a null or a spacelike geodesic, 
respectively, and s^b is the geodesic distance between and x^, i.e. 



SAB = / ^Jlgfj^udxt^dx"] 



The relevance of the world function in the problems related to the light deflection, 
the time delay or the gravitational frequency shift is justified by the following 
properties, which are easily deduced from the above-mentioned theory. 

Property 1. The covariant components of the vectors tangent to the geodesic 
path Tab o,t xa cmd xb respectively, are given by 

9.^^)=-^}-a,xb). (22) 

These fundamental formulae are immediately deduced from equations and 
P5|l . They show that the vectors tangent to the geodesic path V ab at xa and xb can 
be explicitly determined when the world function ^(xat^b) is known. 

Any other affine parameter C along V ab is such that 

c = (Cs - a)A + a , (24) 

where C,a and C,b are the values corresponding to xa and xb^ respectively. As a 
consequence, we have for the tangent vector dx^ jdC, the general formulae 

dx" \ 1 dVL 



9p,u-Tir\ =-7 ^^tt^^a.xb), (25) 

dC J A Cb-Ca ox'X 

dx''\ 1 dn , 

gt^u—rr] =7 — -^{xa,xb)- (26) 

. dC J B Cb-Ca dx% 

It is clear that point xb may be replaced in (|26|l by any point x{Q on Tab which 
differs from xa- Returning to A for the sake of simplicity, we find that the covariant 
components of the vector tangent to Tab at point x{X) are given by 

dx''\ 1 dQ , 

where dVl/dx^ denotes the partial derivative of {1{xatx) with respect to at point 
X. It immediately follows from equations (|27|l that the system of equations 

can be regarded as the first-order differential system governing the geodesic paths 
passing through a given point xa- The regularity of this system at xa is a direct 
consequence of the following property [19]. 

Property 2. The first- order partial derivatives of VL{xa,x) may he expanded as 
— (xA,a;) = g^,i,{xA){x'' - x^a) + Cf,af3{.XA,x){x°' - x'X){x^ - x^) , (29) 
where the functions C^ap{xA,x) remain hounded in the neighbourhood of xa- 
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It results from this property that dfl /dx^ {xa ,x)^0 and that the 

r.h.s. of equations H27|) and (|28|l remains bounded as A ^ 0. 

Now, the foUowing statement can be straightforwardly derived from equations 
iP|l . (O and (|T7|). 

Property 3. The world function VI{xa,xb) satisfies the Hamilton- J acobi 
equations 

^g"^{xA)-^{xA,XB)-^{xA,XB) = QixA^Xs) , (30) 

7:9°''^ixB)-^{xA,XB)-^{xA,XB) = n{xA,XB) ■ (31) 
^ OXg dx'g 

As we shall see below, these equations and Property 2 enable to construct the 
world function in any post-Minkowskian approximation. 

Properties 1, 2 and 3 are valid whatever the nature of the geodesic curve joining 
Xa ant xb- In the case of null geodesies, (|2()|l and L = Ha = Hb = immediately 
lead to the following statement. 

Property 4. Two points xa and xb are joined by a light ray if and only if the 
condition 

nixA,XB)^0 (32) 

is fulfilled. 

Thus, fl{xA, x) — is the equation of the light cone Cx^- 



3. Time transfer functions 



It follows from Property 4 that ifQ{xA,XB) is known in Mxa ^ it is possible to determine 
the travel time tB — tA of a photon connecting two points xa and xb as a function of 
tA, Xa and xb or as a function of t^, xa and xb- It must be pointed out, however, 
that solving the equation fl(ctA, xa, ctB, xb) — for yields two distinct solutions 
and tg since the timelike curve = cuts the light cone Cxy^ at two points x^ 
and x~^, Xg being in the future of x~j^. In the present paper, we always consider xa as 
the point of emission of the photon and as the point of reception, and we focus our 
attention on the determination of t^ — tA (clearly, the determination of tg — tA comes 
within the same methodology). For the sake of brevity, we shall henceforth write Ib 
instead of t^. 

In general, tB—tA may be considered either as a function of the instant of emission 
tA, and of Xa, xb, or as a function of the instant of reception tB and of xa, xb- So 
we are led to introduce two distinct (coordinate) time transfer functions 7^ and % 
respectively defined by 

tB - tA ^ %itA,XA,XB) , (33) 

and 

tB -tA=%{tB,XA,XB) ■ (34) 

We shall call %{tA,XA,XB) the emission time transfer function and %-{tB,XA,XB) 
the reception time transfer function. 
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There exist direct relations between the time transfer functions and the 
components of the vector tangent to a null geodesic. Indeed, it results from equations 
H32|l and that whatever Xg, x''^ and x^, one has the relation 



Differentiating this identity with respect to x\ 
seen that the relations 



cTe{tA, Xa, Xb),Xb) = . (35) 

\ and x^, respectively, it is easily 



1 



d% 
dtA 



{tA,XA,XB) 



= 0, 



dx 



{xa,xb) + c 



A 



on d% 



xa,xb) = 0, 



(36) 



(37) 



dx' 



^{xa,Xb) -g-ritA,XA,XB) + T^ixA,XB) = 



B 



dx\ 



(38) 
Of course, 



hold for any couple of points {xa,xb) connected by a null geodesic 
analogous relations may be derived from the identity 

Vl {x% - cTritB, XA, Xb),XA, x"b,Xb) = . (39) 

Comparing these relations with equations (|22l) - (|23|l . we get the following theorem for 
the components of the vectors tangent to a light ray. 

Theorem 1. Consider a photon emitted at point xa — {ctA,XA) and received at 
point xb — {ctB,XB)- Denote by the vector dx^ /dC, tangent to the photon path, 
C, being any affine parameter. Then, one has relations as follow for the covariant 
components of the vector tangent at xa and x b 



dx\ 



dx 



dx\ 



(fco). 



1 



d% 
dtA 



dtA 



= 1 - 



B 

-1 



1 



d% 

' dt^ 

d%_ 

dx\ 



d% 
dtZ 



(40) 



(41) 



(42) 



where Tp andTj. are taken at (i^i, a;^, a;^) o-nd (tB,XA,XB), respectively. 

These fundamental formulae show that all the theoretical problems related to the 
directions of light rays or to the frequency shifts may be solved as soon as at least 
one of the time transfer functions is explicitly determined. This property will be very 
useful in practice since extracting the time transfer formulae (|33|) or (|34|l from equation 
H32|l . next using equations (|40|) - (|42|l will be more straightforward than deriving the 
vectors tangent at xa and xb from equations H22() - (|23|) . next imposing constraint H32() . 

Case of a stationary space-time - If space-time is stationary, we can choose 
coordinates {x^) such that the metric does not depend on a;". Then, the world function 
is a function of x^ — x^, xa and xb, and the two time transfer functions previously 
introduced are identical. So equations H33|l and H34|) reduce to a single relation of the 
form 



tB -tA = T{XA,XI 



(43) 
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Since T only depends on Xj^ and Xg, it is immediately deduced from equations 
H40|l - (|42|) that the vectors {k^)A and (k^)B defined by their covariant components 

~ dT 

(fco)^ = 1, {ki)A^ c——{xa,xb) , (44) 
~ ~ dT 

(fco)s = 1, {h)B ^ - c—-{xa,xb) (45) 

oxg 

are tangent to any light ray connecting xa and x b , respectively. It must be pointed 
out that these tangent vectors correspond to an affine parameter such that fco = 1 
along the ray (note that such a parameter does not coincide with A). 

4. Determination of VI{xa,xb) in the nth post-Minkowskian approximation 

Equations H31|) or H30|l are the basic formulae enabling to determine the perturbative 
expansion ^ . One might also directly use the integral form of the world function given 
by (|20|l . but the calculations would be very tedious (a brief outline of this method is 
given in Appendix A for the second post-Minkowskian approximation). 

First of all, let us briefly consider Minkowski space-time. The geodesic path 
connecting xa and is curve having parametric equations = a;"p^(A), with 

x'^^^{\)^{x%-x'X)\ + x'\, Q<1<1. (46) 

According to (|46|l . equation (|2U|I gives immediately for the world function in Minkowski 
space-time 

f7(«)(x^, xs) = \vM - x'X){x-'s - x--^) . (47) 

In what follows we suppose that space-time is a perturbation of Minkowski space- 
time. Let us assume for a moment that metric is written in the form 

gfiiy = + . (48) 
Henceforth we suppose that xb € Mx^ ■ The geodesic curve Tab connecting xa 
and is then a perturbation of the straight line F^^. So, the parametric equations 
of Tab may be written in the form 

x°'{l)^{x%-x'X)X + x'X + X"{l), 0<l<l, (49) 
where functions X"' satisfy the boundary conditions 

X"(0) = 0, X"(1) = 0. (50) 
Noting that i'' (Z) = x^ - x^ + X^'H) and that 



V^A^^B-^A)X"il)dl^O 
by virtue of equations ljA.3|) . it may be seen that equation H20() transforms into 

fl{xA,XB) ^ i^''°\xA,XB) + -{^B ^ ^A)i^B ^ ^a) / h^^{x{\))dl 



2 

1 







2 



j^^,X^{l)X'^{l) + 2{x% ~ x^) V(^(A))^"(/) 

+ h^,{x{\))X^'{l)X''{l)\dl. (51) 
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All the integrals in (|51|l are taken over geodesic path Tj^b, which is generally 
unknown. However, it is possible to obtain a perturbative expansion of fl^XA^ xb) 
involving only line integrals taken over the straight line when the metric is given 
by expansion JQ). Indeed, noting that the concomitant expansion of the contravariant 
components (7^"^ is given by 



n=l 



where the set of quantities g^^^ is determined by 

i/(l) — '/ '/ ypa ' 



9(n) 



ri-1 



n—p) ' 



(52) 

(53) 
(54) 



we can state the theorem below. 



Theorem 2. Assuming that the metric is represented by expansion and that 
Xb (zAfxA, the world function is given by equation (|2l, namely 

00 

n{xA, xb) = ixA,XB) + J2 G"^'"^ (^-4, xb) , 



n=l 



where 



1 

n^^\xA,XB) = --^Vt^pVi^aix'^B - x''j^){x''b - x'^a) J 9(i){xio)W)dX, 
n^'''\xA,XB) = -^Vf^pVi^aix'^B - x''a){x% - x'^) j (a:(o) (A))dA 



(55) 



1 

2 Jo 



2 

1 n— 1 

H 9'^p)i^{o)W) 



1 



(Xb-^a)j g^^ ixA,X(^o)W) 

dX 



+ ^^(^^' Ho)W) Q^. i^A, a;(o)(A)) 



9=1 



1 1 dfl^P^ d^^"-~P^ 
^ A2'^'"' ^;^(^'4,2;(o)(A)) (xA,X(o)(A))dA, 

(0) 



(56) 



whatever n > 2. j4// integrals are calculated along straight line defined by equation 

Before we prove this theorem, we can enunciate a corollary as follows. 

Corollary. In terms of covariant components of the metric tensor, Q^^\xa, xb) 
and {xa, Xb) may be written in the form 



n^'\xA,XB) = lix'i, - x'X){x-'b - x-'a) /'gW(x(o)(A))dA , 



n^'H^A,XB) = ^{x^b-A){^b-^a) / 5^,1^ (2^(0) (A)) 



(57) 



-^''"5^pH^(o)(A))5lin^(o)(A)) 
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/■I 1 r)O(l) 



/q a (9a; 
1 1 „„ 



A^''"" ^;^(^'4,a;(o)(A)) ^^(x^, a;(o)(A))dA 



(58) 



This corollary is an immediate consequence of equations IjSSI) and (|56|) for rt = 2 
when equations (|53|I - H54I) are taken into account. 

Proof of Theorem 2. Points xa and xb being given, we may consider n(xA,x) 
as a function of x denoted by Cl{x). This function satisfies the Hamilton- Jacobi 
equation (see equation (p?T|l l 



1 „ dfl d^l — 
_2^"'(^)^(^)a^(^)-"(^)- 
Function may be represented by a series 

oo 

n(x) = n^°\xA,x) + G"l7^"^(x) 



(59) 



(60) 



where each term Q satisfies the boundary condition 

n^"\xA) = o. 



(61) 



Moreover, it follows from Property 2 that the first-order derivatives of each ft must 
satisfy the boundary condition 



lim 

X — >XA 



dxf^ 



(x) 



= 0. 



Inserting developments H52I) and (|60|l into equation (|59(l . and then putting 



9m = 



dxt^ 



■{xa, x) = rii^p{x'^ — Xj^) . 



(62) 

(63) 
(64) 



it is easily seen that equation (|59|) splits up into an infinite set of partial differential 
equations as follows 



dxf" 



where 

^^'Hx) = -^q^q,g>l^^{x), 
$(«)(a^ ^ 



(65) 
(66) 



2 



whatever n > 2. 



E 

9=1 



^x^^ dx'^ 



i9a;'^ 



(67) 
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Assume now that x moves along F^^. This means that x varies as a function of 

A according to equation (|46|l . As a consequence the total derivative of fi^""* along T^^b 
is given by 

^ = «~^A)-Q^i^io)W)^ (68) 



which implies that equation H65|l may be written in the form of an ordinary differential 
equation 

A _ (A)) = (A)) . (69) 

It is easy to show that equation H69|l admits one and only one solution satisfying 
boundary conditions H61(l and 1)62(1 . namely 



(x(o)(A)) = A^ ijCi>(")(a;(o)(A'))dA', (70) 

from which one can deduce the formulae H55|) - (|56|l . The only difficulty is to prove that 
A~^$^"'' (x(o) (A)) is a continuous function of A in the range < A < 1. The property 
is obviously true for n = 1, since it follows from equation (|66|l that 

$(^)(x(o)(A)) = -^{xb - ^Di^B - xa)VppV^<j9^^i){x{o)W) ■ (71) 

As a consequence, fi*'^^ (a;(o) (A)) is given by ((70|l for n ^ 1. Therefore, the formula 
H55|l holds. Replacing a;^ by a; in ((55|l and then differentiating with respect to x gives 







— {xA,x)^-qpj^ ?7^<xff('i'')(a:(o)(A))dA (72) 



dx 

1 







where a;(o)(A) is defined by 

x^o^{\)^{x" -x1)X + x1. (73) 
Setting X = a;(o) (A) in ifT^ yields 

a _ - - 

-^-^{xA:X(^a)W) = -Kx% - X'X) J r^apriuag1l){x(Q){\))d\ 

- iA2(x^ -a;S)(4 -4)^\„,,7^,g[';') ,^(af(o)(A))AdA. (74) 

Inserting this expression into equation 1(67(1 written for n = 2, it is easily seen that 
;v-2(j)(2) jg continuous function of A. As a consequence, fl is given by ((70(1 for 
n = 2, which proves that the formula 1(56(1 is valid for n = 2. The above reasoning 
may be used to prove that A"^^'^"^ is continuous whatever n. Q. E. D. 

The expression of nf^^^XA^XB) given by the formula ((58(1 leads to heavy 
calculations owing to the presence of quadratic terms with respect to the derivatives 
of ri^^-' in one of the integrals. Fortunately, a drastic simplification is possible. 
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Theorem 3. The second-order term fl'^'^^ {xa, xb) may be written in the form 



1 ^ 



gj-H^wW) - ^''^9l';(^(o)(A)).g«(x(o)(A)) 
«q(i) 

-V'^ i^A , a:(o) (A)) gl% (x^o) (A)) 



dX 



1 arid). arid) 

.(0) 



(75) 



the integral being taken along ^aq- 

Proof of Theorem 3. It follows from equation H74() that 

A (^^' ^(0) (a^^' 2:(o) (A)) 



lim 

A-.0 



= 0. 



(76) 



Taking this equation into account, an integration by parts gives 



1 r' 1 



ai7d 



9»d) 

(2:^,2:5) (2:^,2:3) 



;^2'?'"' ^^(2^^'2;(o)(A)) ^^(a;^,a;(o)(A))dA 



2 9a;'; 



dX. 



(77) 



In order to transform the r.h.s. of equation (|77|l . let us perform the change of variable 
A' = A A in equation (|74|l . We obtain 



917( 



9x 



-(xA,X(o)(A)) = (a;^-a;^) / (^.^^^ (A'))dA' 



+ 2(^b-2;^)(^b-^a) / 5^y,.(2;(o)(A'))A'dA'. 



(78) 



Differentiating (|78ll with respect to A yields 



d 

dX 



(9i7d) 

^^(xA,a;(o)(A)) 



(x^ ~x'\)gfJ{x^Q){X)) 



+ 2(^B - - 2^A)A5^'j,.(x(o)(A)) 



(79) 



Substituting this expression into (|77|) , and then inserting the result of this substitution 
into (|55|l . we finally get the formula (|75|l . Q. E. D. 

5. Time transfer functions in the nth post-Minkowskian approximation 

It follows from equations H35|l and (|47|) that function Te{tA,Xa,XB) satisfies the 
equation 



c'T^^ ^RIb-^Y. G''^^''\^a,xa, x\ + c%,xb) , 



(80) 
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where 

Rab = \xb - xa\ ■ 

Since xa is the point of emission and xb the point of reception, we are 
looking for the solution to equation l|Sn|l such that cTe{tA,XA,XB) = Rab + 0{G). 
Mathematically, our problem consists in finding a function u{x,e) governed by the 
equation 

oo 

u2 ^a;2 + ^e"7i„(u) (81) 

and satisfying the condition 

u{x,e)\^^o = X , (82) 

e being a parameter sufficiently small to warrant the convergence of the series which 
are handled. Using Maclaurin's expansion of u about e — 0, the solution to equation 
H81|l may be written in the form 

oo 

u{x, e) — x + e^Un{x) , 

n=l 

where 

1 d^u{x,e) 



Unix) 



nl de" 



e=0 

Performing successive differentiations of (|81|l with respect to e and then taking (|82|l 
into account, a straightforward calculation yields Un{x) for any n. The first terms are 
given by 

hi{x) 

U2{x) = — [h2ix) + h[{x)ui{x) - ul{x)] , 



2x 



h^{x) + h'2{x)ui{x) + -h'l{x)u\{x) + h\{x)u2{x) — 2ui{x)u2{x) 



(83) 

Applying this procedure to equation (jHUJ, it would be possible to get 7^ at any 
post-Minkowskian order. For the sake of brevity, let us define the following functions 

of tA,XA,XB: 

ni''HtA,XA,XB)^n'^"HA.XA,x''A + RAB,XB) (84) 

and 

~(n) cj'^rj^") 

^e\k(*A,XA,XB) = , {x°^,Xa,X°a + RaB,Xb) (85) 

for k = 1,2,3, ... (be careful that is not a partial derivative of fli"'^). Applying 

equations H83() . we can state the following theorem. 

Theorem 4. In the third post-Minkowskian approximation the emission time 
transfer function is given by 
1 3 

%itA,XA,XB) = -RAB + Y.'^"'^e"^(*A,XA,XB) + 0(G^), (86) 
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(1) 



r(2) 



cRab 
1 

cRab 
1 

cRab 
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(87) 



(89) 

(90) 
(91) 



the functions fii"^ and ^^J^^j^ being defined by (|84() and (|85|l , respectively. 
Of course, defining 

f2(") (ts , a;^, a^s) = (x^ - , a^A, , a;^) 

and 

~(n) 

fi^|fc(ii3,a;A,a;s) = (ga;0 )fc (^a ^ -Rab, a^A, a;^, aJs) , 
a similar theorem can be stated for Tj-ltg, xa, xb)- 

Theorem 4 a. In the third post-Minkowskian approximation the reception time 
transfer function is given by 

1 ^ 

Tr{tB,XA,XB) = -RaB + y^G''r}^\tB.XA.XB) + 0{G^) , (92) 

where 



n 



(1) 



cRai 



r(^) = -- 



cRab 



cRab 



(93) 
(94) 

(95) 

It may be seen that Tr^\ 7^'-"'' can be determined once Te'^\ T^*-"-* are 
known and conversely. Let us examine in detail the cases where n = I and n = 2. 
Comparing expressions H87|l and (|93|l . we get immediately 



and conversely 



X^'\t,XA.XB)^Tj''> (^t-^,XA,XB^ 

rW(i,a;A,a=s)=rW (^t+^,XA,XB^ 



(96) 



(97) 



Now comparing expressions H88|) and H94|l . and then noting that (|87|l implies the 
identity 

OTP 



dt 



-{t,XA,XB) 
1 



Rab 



{x^,XA,X° + Rab,xb) + ^^^{x°,xa,x" + Rab,xi 



dx\ 



dx% 



(98) 
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we get 



Rai 



%^^\tB^XA,XB) = TP [t 

Conversely, we have 

= r}') [tA 



Rai 



-,xa,xb 



-,Xa,Xb 



dTe 



(1) 



dt 



Rai 



,Xa,Xb 



, ^Rab 

lA H ,Xa,Xb 



dV 



(1) 



dt 



Rab 

ts ,Xa,Xb 



. ^Rab 

tA H ,Xa,Xb 
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(99) 



(100) 



6. Application to a static, spherically symmetric body 

In order to apply the general results obtained in sections 4 and 5, let us determine 

the world function and the time transfer function in the second post-Minkowskian 
approximation. Choosing spatial isotropic coordinates, we suppose that the metric 
components may be written as 



(2) 

5oo 



2M 



2(3- 



(2) 
90^ 



0, 



9 



(101) 



0, 



(2) 3 M2 



2 c^r^ ' 

where /3 and 7 are the usual post-Newtonian parameters and 5 is a post-post- 
Newtonian parameter (in general relativity, /3 = 7 = (5 = 1). Furthermore, we suppose 
that XA and xb are such that the connecting geodesic path is both entirely outside 
the body and in the weak-field region. We use the notations 

r = \x\, rA = \xA\, tb^IxbI- 

It is easily seen that GQ'^^\xa,xb) is the term due to the mass in the multipole 
expansion of the world function given by equation (56) in [22] (see also [31]) 

ft^'Kx'i, XA, x%,xb) = -^ [ix% - x^Af + 7-Rab] FixA, xb) , 
where 



(102) 



F{xa,xb) 



dX 



|a;(o)(A)| Rab 



In 



For the integrals involving gj^J and terms quadratic in gi^V in the expression of 



TA + T-B + Rab 
ta+tb- Rab 
(1) 



(103) 



Q.^'^\xa,xb), we find 
-(a;^ — x^){xb — Xa) 

= ^[(/3~2)(x°,-x^)^+(27^-^)i?: 



5^in^(0)(A)) - r?^'^5iy(^(o)(A))5«(a;(o)(A)) 



dX 



ab 



E{xa,xb) , 



(104) 



where E{xa,xb) is defined as 
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1 



— arctan 
arccos(n^ • ub) 



{xb-xa)-xb 
arctan — 

VM^I - i^A ■ XbY 
{xb - xa) ■ XA 



V^a'^b ~ i^A ■ xb) 



rArB\Jl - [riA ■ nB^ 
riA and Ub being defined as 

UA = XA/rA , TlB = Xb/tb ■ 

For the integral involving the gradient of fi^^^ in the formula we obtain 



(105) 



-{x%~x^^){xI~x'^a) / V'" 



dxP 



{XA,X(Q) (A)) (^x^Q) {\))d\ 



1 + ^ 



1 {X% - + -fH 



AB 



-2-fE{XA,XB) 



TB J Rab '^atb + {xA ■ Xb) 
{xl-x\f -^R\sF{xa.xb) 



R-AB 



TB 



(106) 

Substituting p04|) and pOt)|) into H75|) . and then carrying out the calculation of 
the square of the gradient of il^^^XA, xb) at xb, we obtain an expression for the world 
function as follows 



nixA,XB)^^ix%-x"^f -^R\s 



GM 



[{X%~x\f+^R\B\F{XA.XB) 



GHP ( {x' 



^AB 



rATB + [xA ■ Xb) 2 



]-F'^{xa,xb) 



27 



tatb + {xA ■ Xi 



- {2- /3 + 2-f)E{xA,XB) 



+ i2 + -f)F\xA,XB) 



+R 



AB 



7^ 3(5 7^ 

\ ~ -rF{xA,XB) - —F^{xa,xb) 

tatb + [xA ■ Xb) 4 2 



+ 0{G^) . 



(107) 

Now, using the formulae (|87|l and l|88|l . a straightforward calculation gives for the 
time transfer function 

^, , Rab , (7 + l)G'A/ frA+rB+RAB 
T{xa,xb) ^ H ^ In 



G^M^Rab 



\rA + rB - Rai 
4/3 + 87 + 3(5) arccos(n^ ■ tib) 
W'^A^B - i^A ■ xb)^ 



(1 + 7)' 



rA tb + {xA ■ Xb) 



0{G^ 



(108) 

The first-order term is the Shapiro time delay [37]. The second-order terms in 
H108fl generalize the result found by Brumberg in the case of general relativity [32]. 
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Denoting by ^ab the angle formed by vectors and and noticing that 
E{xa,Xb) reads 

arccos(n^ ■ ub) 4'ab 



■s/r\ r| - {XA ■ XbY rArBsinipAB 
we find an expression for T(xa,xb) as follows 



< tpAB < TT , 



T{xa,xb) 



Rab 



(7 + l)G M frA +rB + RA 



■In 



TA 



c^rATB 



(2-/3 + 27+^5).^ 
4 'sva^AB 



rB - Rab 



1 + cos %1}AB 

which may be very helpful for numerical estimates. 



(109) 



0(G3), 



(110) 



7. Conclusion 



The central result of this paper is given by Theorem 2. Assuming that xb is in 
the normal convex neighbourhood of xa^ we have shown that U,{xa,xb) can be 
obtained in the nth post-Minkowskian approximation by a recursive procedure which 
spares the trouble of solving the geodesic equations. Any nth-order perturbation 
^'''^\xa,xb) is an integral taken along the zeroth-order straight line joining xa 
and Xb- Moreover, in Theorem 3 we have found a remarkable simplification of the 
integral giving n'^'^\xA,XB)- In Theor ems 4 and 4a we have outlined a recursive 
procedure enabling to obtain the perturbations of the time transfer functions when the 
perturbations of U,{xaiXb) are known, assuming that the null cone at xa is included 
in the normal convex neighbourhood of xa- Since the time transfer functions are 
sufficient to determine the direction of a light ray in the 3-space relative to a given 
observer (see Theorem 1), the systematic methods explored here seem very promising 
to tackle the relativistic problems raised by highly accurate astrometry in space. 

Using the simplification found in the second post-Minkowskian approximation, 
we have shown that the calculation of U!^"^^ {xa,xb) and T'^'^\xatXb) reduces to 
elementary integrations in the case of a static, spherically symmetric space-time. The 
simplification within the second order works out so well that the question whether 
it works also for the higher-order terms naturally arises. We hope to return to this 
problem later. For the time being, we are applying the results obtained in the static, 
spherically symmetric case to the description of the light deflection between two points 
located at a finite distance. 



Appendix A. Another determination of n'^'^^ {xa,xb) 

It is possible to derive ^^'^\xatXb) directly from the Euler-Lagrange equations of 
geodesies, namely 

d (, , ,dx''\ 1, dxi" dx" , 
^(^(.W + V)^j - 2'^.-.^^- (A.l) 

It is clear that the functions AT'' defined by equations H49() admit the following 
expansions 

oc 

X^(A,G) = ^G"A[;)(A), (A.2) 

n=l 
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where the functions ^[^-j W satisfy the boundary conditions 

X^^^iO)^0, X^'^^W^O Vn>l. (A.3) 

Substituting for gf^^, and from equations and l|A.2|l respectively into 
equation (|5T|l gives 

1 



1 ^ 



Noting that 

g^^^Ml)) - 3l'.^(^(o)(A)) + GX(;)gW^(x(o)(A)) + 0(0^) 
wc find that il*^^^ and fi'^^ are given by 

^'^'\xa.Xb) = Wb-X^j,){xI-x\) [\l^J{xrM)dX, 



1 



(A.4) 
(A.5) 

(A.6) 



+2(a;^ - x:;)i:fi)(Og(i,n^(o)(A)) + 7?^.i:(;)(0i:(;)(0} rfA , 



(A.7) 



where the integrals are now calculated along unperturbed geodesic F^^. 

Equation ljA.6|) coincides with the formula (|57|l . Let us now transform the r.h.s. 
of equation (|A.7p in order to get rid of all the terms involving X'^^-^ or X'^^y Using 
equations \A.\\ and retaining only the terms of first order in G, we obtain 

~ x^^){x'^B - x''A)9%{H^m = ^ [^P-^ri)(0 + [x-h - x^)g(;.n^(o)(A))' 

Substituting the r.h.s. of these equations into equation (jA.7|l . and then integrating by 
parts, we find 



n^^\xA.XB) = \l {x^B- Oi^B - 2^^)511^(2^(0) (A)) - r^,.X'^,^{l)Xl,^{l) 



dX. 



In order to transform the integral 



1 



^(2) 



rj^,X^^^{l)X^,^{l)dX, 



(A.. 



(A.9) 



let us insert expansions (|A.2p and equation Q into equations (I28|l . Retaining only 
the terms of first order in G, we find that A^^-i satisfies the differential equation 



^(1) " A^(i) - A 



1 r}n(i) 

v'"'^{xa,x^o){1))~v'''{x^b~^a)9^^Xho)W) 



(A.IO) 
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dX' . 



The only solution to this equation satisfying the boundary condition X^^~^{1) = is 
given by 

X^^^iX) = XYf^^iX) , (A.ll) 

where 

-{xg - XA)g)J{xi^a){X )) - ^-q^^ 

(A.12) 

Using Property 2, it is easily checked that this solution also satisfies the boundary 
condition ^[^^(O) = 0. 

Now substituting for X^^j^^ (A) from X'^^s^ (A) = XY^^^ (A) + y^*^-) (A) into equation 
"9|) and then integrating by parts, we get 



^(2) 



X^il^.Yl{^{\)Y^,){X)dX., 



(A.13) 



since Y^^-^{1) — 0. Differentiating the r.h.s. of equations ljA.12|) with respect to A and 
then substituting into equation (|A.13p . we finally obtain 

5(2) - \{^% - A){^B - f r?'"^5«(^(o)(A)).g«(a;(o)(A))dA 



' 1 ..M^'> 



O TV"" 



A' dx" 



{xA,X(Q) (A))^^ (^x^^^) {X))dX 



(xA,a;(o)(A)) -—-;^{xA,x^^Q){X))dX. 



2 Jo A2 ' dxt" ^"--'-^ww/ 
Substituting for S(2) from this last equation into equation (jA.8|l . we recover the 
expression of fl'-^^ {xa, xb) given by the formula (|58|l . 

As a final remark, we shall note that equations l|A.ll|l and l|A.12|l yield the first- 
order perturbative term of the geodesic path joining two given points xa and xb as 
integrals taken over the straight line with ends at xa and xb- 
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